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Abstract 

In this paper we discuss the question of how to decide when a general chemical 
reaction system is incapable of admitting multiple equilibria, regardless of param- 
eter values such as reaction rate constants, and regardless of the type of chemical 
kinetics, such as mass-action kinetics, Michaelis-Menten kinetics, etc. Our results 
relate previously described linear algebraic and graph-theoretic conditions for in- 
jectivity of chemical reaction systems. After developing a translation between the 
two formalisms, we show that a graph-theoretic test developed earlier in the con- 
text of systems with mass action kinetics, can be applied to reaction systems with 
arbitrary kinetics. The test, which is easy to implement algorithmically, and can 
often be decided without the need for any computation, rules out the possibility of 
multiple equilibria for the systems in question. 
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1 Introductory material 



There is increasing interest in methodologies for drawing conclusions about 
the dynamics of a chemical reaction network based only on the network struc- 
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ture, i.e., with limited or absent knowledge of the kinetics. Early work in this 
direction is exemplified by [T3lll0fl2|[TT] . with more recent strands including 
discussions of monotonicity P^IPIT] . and discussions of injectivity [5117115] . Al- 
though the reaction systems discussed in examples are often assumed to have 
mass action kinetics, an important feature of PIT|3] was that they described 
criteria which applied for more general chemical kinetics. The aim of this paper 
is to extend graph-theoretic results, which were developed in [7] in the context 
of mass action kinetics, to the case of general chemical kinetics. These graph- 
theoretic criteria are more restrictive than the matrix-theoretic results in [3J, 
but are more intuitive, and give rise to conditions which are less expensive 
computationally, and are often easy to check by hand. 

Dynamical systems derived from chemical reaction networks. A chem- 
ical reaction system in which n reactants participate in m reactions has dy- 
namics governed by the ordinary differential equation 

X = Sv{x), (1) 

where x = [xi, . . . ,Xn\'^ is the nonnegative n- vector of reactant concentra- 
tions, V = [f 1, . . . , Vm]'^ is the m-vector of reaction rates and S is the n x m 
stoichiometric matrix. Arbitrary orderings can be chosen on the sets of sub- 
strates and reactions. Further, 5* is only defined up to an arbitrary re-signing 
of its columns, equivalent to a switching of the left and right hand sides of a 
reaction. It is trivial that all results here are independent of the orders chosen 
on substrates and reactions. We will also confirm below that they are indepen- 
dent of the signing of columns of S. System ([1]) defines a dynamical system on 
the nonnegative orthant in M". With the additional assumption that all sub- 
stances may have some infiow (which is allowed to be zero) and some outfiow 
which increases with concentration, we obtain the related system 

X = K + Sv{x) -Q{x). (2) 

Here i^' is a nonnegative vector representing the infiows and the diagonal 
function Q{x) = [Qi{xi), . . . , Qnixn)]^ represents the outfiows, and is assumed 
to satisfy > for each i. The system has Jacobian J = SV{x) — D{x) 
where the m x n matrix V{x) is defined by Vij{x) = and the diagonal 

matrix D{x) is defined by Dai^x) = It is notationally convenient to omit 
the exphcit dependence on x, and write J = SV — D. 

Previous results on injectivity and the main result of this paper. In |3] 

and [1] a reaction system was termed "nonautocatalytic" (NAC for short) if 5* 
and V'^ have opposite sign structures in the following sense: SijVji < for all i 
and j, and Sij = ^ V,j = 0. These conditions are naturally fulfilled provided 
that no substrate occurs on both sides of a reaction (either with the same or 
with a changed stoichiometry) . Here we refer to systems where the above two 
conditions are fulfilled, as NIC reaction systems, in order to emphasise that 
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the conditions only rule out one-step catalysis. Note that the NIC condition 
is not very restrictive for realistic biochemical reaction networks: for example, 
all the enzymatic reactions considered in [8| satisfy this condition. 

A square matrix M is sign nonsingular if all matrices with the same sign pat- 
tern as M are nonsingular [4j. In [3j it was shown that System ([2]) is injective, 
and hence incapable of multiple equilibria, provided that the stoichiometric 
matrix S is strongly sign determined (SSD), i.e., all square submatrices 
of S are either sign nonsingular or singular. This led to a characterization of 
injectivity based entirely on a computation on the stoichiometric matrix. On 
the other hand, following theory developed in [5], a signed, labelled, bipartite 
multigraph termed the Species-Reaction graph or SR graph was constructed 
in and used to make claims about the nonexistence of multiple equilibria 
in systems of chemical reactions with mass-action kinetics. This time, rather 
than a matrix computation, a certain graph-theoretic condition, sometimes 
checkable by observation alone without the need for any computation, was 
shown to be sufficient to guarantee the absence of multiple equilibria. Here we 
will combine these ideas, and in fact show that the graph-theoretic condition 
in [7] suffices to guarantee the absence of multiple equilibria for NIC systems 
with arbitrary kinetics. 

We will define an SR graph for an arbitrary NIC reaction system solely via 
the stoichiometric matrix 5* of the system. We will then present a condition 
on this graph which will be termed Condition (*). Our main result is the 
following: 

Theorem 1 Consider the SR graph G of an NIC reaction system with stoi- 
chiometric matrix S. Assume that Condition (*) is fulfilled. Then S is an SSD 
matrix. 

As a consequence, if Condition (*) holds for the SR graph associated with 
an NIC reaction system, then System ([2]), i.e. the system with outflows, is 
injective, and hence incapable of multiple equilibria. Along the way to proving 
Theorem [T] we will prove several auxiliary results relating computations on 
matrices to computations on graphs. Of these. Corollary [13] is of independent 
interest . 



2 The SR graph 

The SR graph, introduced in [7j, is a bipartite graph with n substrate vertices 
or S-vertices, m reaction vertices or R-vertices, and an edge between S- 
vertex i and R- vertex j iff substrate i participates in reaction j. For an NIC 
reaction system there is a one-to-one correspondence between edges in the SR 
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graph and nonzero entries in the stoichiometric matrix S: clearly Sij im- 
plies that substrate i participates in reaction j, and given the NIC condition 
it can participate on one side of the reaction only, giving one edge between 
S- vertex i and R- vertex j (in the general case multiple edges are allowed). 
On the other hand, given the NIC condition Sij = implies that substrate i 
does not participate in reaction j, and hence that there is no edge connecting 
S- vertex i and R- vertex j. 

In this paper, for convenience, we make slight modifications to the labelling 
of the SR graph in [7j. An edge Cij (where i E {1, . . . ,n} and j G {1, . . . , m}) 
will be taken to mean an edge between the ith S-vertex and the jth R- vertex. 
Such an edge exists if and only if Sij ^ 0. Rather than labelling edge Cij with 
the complex label associated with substrate i and reaction j as done in [7j, 
we simply give it a sign, so that sign(ejj) = sign(5'ij). Since S is only defined 
upto an arbitrary signing its columns, the signing of edges in the SR graph 
is not unique. However, as we will show in Lemma [2] after developing some 
further ideas, all results here are independent of the choice of signing. Whether 
labelled with the complex labels or signed, the important thing is that the SR 
graph contains information on whether substrates occur on the same side of 
a reaction or on opposite sides of a reaction. 

Given the one-to-one correspondence 

Sij *v=7^ Cij 

when Sij 7^ it is convenient to allow a slight abuse of notation and refer to 
Sij as an edge in G. Two edges S'j^j^ and Si^j^ share an S-vertex when ii = i2, 
and they share an R- vertex when ji = j2 . 

Remark. As a consequence of our formal redefinition, rather than being asso- 
ciated directly with chemical reaction networks, SR graphs are now associated 
directly with real matrices. 

Example. An example of a matrix and its corresponding SR graph is shown 
in Figure m 
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Fig. 1. Left. A matrix M. Assume that a,b,c,d,e, f,g,h,j > 0. Right. The corre- 
sponding SR graph. S-vertices have been labelled 81,82,33 while R- vertices have 
been labelled Ri, R2, R3, Ra. Positive edges are bold lines while negative edges are 
dashed lines. These conventions will be followed in drawings of SR graphs. 
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Terminology and definitions. Given a graph G, we will say that some ver- 
tex/edge lies in G meaning that it lies in the vertex/edge-set of G. We will 
generally define subgraphs by their edge-sets, i.e., writing E = {ei, . . . ,6^} 
where Cj are edges in G, will mean that E is the subgraph of G containing 
exactly edges ei, . . . , and the vertices on which they are incident. The terms 
path and cycle will be taken to mean open/closed simple walks in G respec- 
tively and are particular examples of subgraphs of G. In the usual way, the 
size of a subgraph E is the number of edges in E, written \E\. When E is a 
cycle or a path this will also be called the length of E. Because of the bipartite 
nature of SR graphs, all cycles are of even length including alternate S- and 
R-vertices. If two paths or cycles are edge-disjoint they may still share some 
vertices. When they share no vertices they will be termed vertex-disjoint. A 
path between an S-vertex and an R-vertex is called an S-to-R path. Note 
that S-to-R paths are of odd length. 

Cycles in an SR graph have a natural parity - they are either odd or even. 
We use a definition equivalent to that in but using the labelling defined 
above. 

Consider a matrix S and the corresponding SR graph G. Let E be any sub- 
graph of G. The sign of an edge e in G has already been defined above, from 
which we can define the sign of E to be 

sign(£') = Yl sign (e) . 

When is even, we define the parity of E to be 

P(E) = (-l)l^l/2sign(E). 

Since cycles are always of even length, the parity of a cycle is always defined. A 
cycle C will be termed an e-cycle if P{C) = 1 and an o-cycle if P{C) = — 1. 
Note that by these definitions, for an e-cycle C we have 

sign(C) = 

and similarly for an o-cycle C we have 

sign(C) = (-1)1^1/2-1. 

We define the value of edge Cij to be val(ejj) = \Sij\, and for a subgraph E, 

val(E) = n val(e) 

When C is a cycle containing edges Ci, 62, . . . , such that and Cj+i mod 2r 
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are adjacent for each i = 1, 



, 2r, we can define 



stoich(C) 



i=l 



i=l 



Note that this definition is independent of the starting point chosen on the 
cycle. A cycle with stoich(C) = is termed an s-cycle. This definition is 
equivalent to the definition in [7]. 

The intersection of two cycles in an SR graph can be divided into a set of 
vertex-disjoint paths. We say that two cycles have S-to-R intersection, if 
each component of their intersection is an S-to-R path, i.e. a path between an 
S-vertex and an R-vertex. 



We finish this section with a lemma confirming that re-signing the columns of 
the stoichiometric matrix does not alter the nature of cycles in an SR graph. 

Lemma 2 Consider a matrix S and the associated SR graph G. Now consider 
any re-signed version of S, say S with corresponding SR graph G . Cycles in 
G are in one-to-one correspondence with those in G . e-cycles (o-cycles) in G 
correspond to e-cycles (o-cycles) in G . s-cycles in G correspond to s-cycles in 
G'. 



PROOF. Ignoring the signs on edges, clearly G and G' are identical graphs, 
so cycles in G and G are in one-to-one correspondence. Consider some cycle 
C E G and the corresponding cycle G E G . Re-signing column j of S means 
re-signing all edges incident on R-vertex j in G . But clearly G, and hence G , 
contains either 2 or edges incident on R-vertex j. So the re-signing causes 
no change in sign(C), i.e., sign(C) = sign(C ), so that e-cycles and o-cycles 
are preserved. Finally, re-signing a column does not alter any of the values of 
edges, leaving stoich(C) unchanged. Thus the re-signing does not affect the 
property of a cycle being an s-cycle. □ 



As a final note, it is obvious that S-to-R paths remain so after a re-signing of 
the graph. 



3 Determinants, permutations and cycles 



Consider some n x m matrix 5* and the associated SR graph G. 

Notation. S{'~f\6) will refer to the submatrix of S with rows from some set 
7 C {1, . . . , n} and columns from some set 5 C {1, ... , m}. If S^jIS) is square, 
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then its determinant will be written 5'[7|(5]. Each submatrix S{'~f\S) corresponds 
to a subgraph of G which we will term G{'-f\6). 

Determinant expansions and term subgraphs. Consider any sets 7 = 
{71, 72, ■■■,7k} C {1, . . . , n} and (5 = {61,62, ■■■, 6k} C {1, . • • , m}, choosing 
7i < 72 < • • • < 7fc and 5i < ^2 < • • • < ^fc, so that both sets have a natural 
ordering. Consider the square submatrix 3(^16) of S. Any permutation a of 5 
gives a term in the expansion of 5'[7|5] of the form Tq, = P(a)5'^^Qj ■ ■ ■ S^^aj,, 
where P{a) = — 1 if a is an odd permutation, and P{a) = 1 otherwise. Note 
that Tq, is simply a real number. Thus terms in the determinant expansion of 
S{'y\6) are in one-to-one correspondence with permutations of 6. 

If To, 7^ 0, then a can also be identified with a subgraph of 0(^^16), 

Each S-vertex chosen from 7 and each R-vertex from 6 occurs in exactly one 
edge in Ea- This follows since each member of {71,..., 7^} occurs exactly 
once as a first subscript in S^y-^^ai, ■ ■ ■ , S^^^ak and similarly each of {61, ... ,6^} 
occurs exactly once as a second subscript in this expression. As a result no 
two edges in Ea share a vertex, and Ea can contain no cycles. We will refer to 
a subgraph of G{'^\6) with these properties as a term subgraph of G(7|(5). 
Clearly term subgraphs are in one-to-one correspondence with nonzero terms 
in the expansion of 5'[7|5]. 

Permutations of a fixed set 6 form a group, and so it makes sense to talk 
about the operations of composition and inversion. Given permutations a, (3 
we will often be interested in the permutation (3 o which takes a to (3. 
Permutations can be written as products of disjoint cycles. A nontrivial cycle 
will refer to a cycle of length greater than 1. Below we will show that there 
is a close relationship between cycles in the decomposition of a permutation, 
and cycles in the SR graph. Throughout this paper, the word "cycle" has 
two distinct meanings: either a special kind of subgraph in an SR graph, or 
a special kind of permutation. In general the meaning will be clear from the 
context. 

Several key constructions in this paper rely on taking two term subgraphs 
corresponding to two distinct permutations of 5, say a and /3, and looking 
at the structure of their union Ea U Ejj. Consider this union: the S-vertex 7^ 
occurs in exactly one edge in Ea U Ep if = and exactly two edges in 
Ea U Ep otherwise. On the other hand, assume 6j = = (3k'. if i = k, then the 
R-vertex 6j occurs in exactly one edge in Ea U i?^ while if i ^ k then it occurs 
in exactly two edges in EaUEjj. Thus Ea^Ep consists of a set of components, 
each of which is either i) an isolated edge corresponding to a trivial cycle in 
j3oa~^, or ii) a cycle corresponding to a nontrivial cycle in f3oa~^. The explicit 
construction is carried out in Lemma [3] below. Any given vertex in Ea U Ej^ 
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has exactly one edge from Ea U E/j incident on it if it corresponds to a trivial 
cycle in /5 o a~^, and exactly two edges incident on it if it corresponds to a 
nontrivial cycle. In a similar way, given a set of k permutations ai, . . . ,ak, 
each vertex in U^Lj^ii^a. has between 1 and k edges from u'^^^^^Ea. incident on 
it. 

Notation. Given a set 7 = {71, 72, . . . , 7/,.}, we write 7 = [71, 72, . . . , 7fe], 
when it is important to stress that 7 is an ordered set. When a permutation is 
written as a product of cycles, we use round brackets to denote these cycles. 
These may include or exclude trivial cycles. 

Example. Let 

S ^ [61,62,63,64], a ^ [61,63,64,62], P ^ [62,64,63,61]. 

Written as products of disjoint cycles, a = (5i)(52, ^3, ^4), /3 = ((^3)(^i, ^2, (^4) 
and /3 o q;-^ = {61, 62) {63, 64). Finally P{a) = P{P) = +1 and P(/3 o a'^) = 
P{P)P{a) = +1. 

Notation. Given that cycles are the fundamental objects in the theory being 
developed in this paper, when we have an index k which is known to belong to 
a set {1, . . . , r}, counting is always done on a circle of size r, so that k+p means 
{k + p — 1 mod r) + 1. Adopting this convention avoids lengthy subscripts. 

The next lemma illustrates the relationship between cycles in permutations 
and cycles in SR graphs, and also begins the process of hnking statements 
about SR graphs and stoichiometric matrices. 

Lemma 3 Consider a matrix S and the associated SR graph G. Assume that 
there is some square submatrix S{'y\6) and permutations a and P of 6 such that 

Ta and Tp are nonzero terms in S['j\S]. Then corresponding to each nontrivial 
cycle in the decomposition fi o there is a cycle in G. In particular G 
contains at least one cycle. 



PROOF. The two nonzero terms can be written exphcitly as: 

Ta = P{Q.)S^f^ai ■ ■ ■ S-y^a^, TfS = P{(3)Sy^p^ ■ ■ ■ S^^p^ . 

Write 13 o a^^ as the product of disjoint cycles. Since (3 ^ a this product 
contains at least one nontrivial cycle. Consider such a cycle of length r (2 < 
r < k), (^6(1), . . . , 5fe(r))- For each j G {!,..., r} define the index a(j) by 
(^a{3) — ^b{j)- (Note that since a is a permutation, ji ^ j2 cbiji) 7^ <^{j2)-) 
The existence of the cycle means that for each j e {1, . . . , r}, S^^^^^s^j) £ -E'a 
and -S^^yj^jy+i) e Efj. Writing these as 
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makes it clear that they form a cycle of length 2r in G. 



□ 



The construction in Lemma [3] will be used frequently - i.e., given two different 
permutations a and /? of a set 5, cycles in /? o will be used to infer the 
existence of index sets {a(j)} and {&(j)} and corresponding cycles in the SR 
graph. When using the construction, for notational brevity we will write ai 
for 7a(i) and 61 for 5f,(i). 

Example. Let 

7 = [71, 72, 73, 74] 

^ = [(^1, (^2, ^3, (^4] 

a = [5i, 53, 54, 82] 

(3 = [§2, ^4, ^3, ^l] 

Let 5(715) be a submatrix of a matrix S and G'(7|5) the associated subgraph. 
a and /5 define the term subgraphs of G'(7|5): = {S-y-^Si, S^^Ss^ 'S'7354, S^^S2} 
and i?^ = {S^^S2y S^^Si^ ^"/sSsy S.y^Si}- From the previous example, /? o = 
(f^i, '^2)(53, ^4), so there are two cycles Ci and C2 in U Ef^. Written out 
explicitly, these cycles are 

C*i = {5*7^5^, 5*^152, *S''y452, ^■y4Si} s-^d C2 = {S^^s-ij 3^264^ 'S'7354, ^ri^}- 



So far we have focussed on constructing cycles in an SR graph from pairs 
of nonzero terms in a determinant. However the reverse is also important 
- inferring the existence of pairs of nonzero terms in a determinant from 
structures in a graph. The basic operation which allows us to do this is a 
particular splitting of a cycle. Any cycle 

can be uniquely partitioned into two vertex-disjoint subgraphs of equal size: 

C = {Saibi, Sa2b2: • • • , Sa^h^} U {Saib2: 'S'aafea, • • • ' ^arbi}- (3) 

We will call this a disconnecting partition of C. With this notion, confirm- 
ing if a cycle C is an s-cycle now involves: 

(1) Constructing a disconnecting partition of C into {Ci,C2}, 

(2) Confirming that val(Ci) = val(C2). 
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4 Preliminciry lemmEis 



With the machinery set up above, we are ready to prove some lemmas. The 
idea of these lemmas is that each of them is quite brief and reusable, so that 
the proofs of the main results in this paper become simpler. 

The first lemma gives us a basic way of checking whether a permutation, 
written as a product of cycles, is even or odd. 

Lemma 4 Consider a permutation a written as a product of cycles from some 
set C. Let 9 = Ucecc Then 



i.e., a is an even or odd permutation according to whether the total number 
of elements in cycles, minus the number of cycles is even or odd. 

PROOF. This follows by writing any permutation as the product of disjoint 
cycles and noting the elementary result that a /c-cycle is an even permutation 
if k is odd and vice versa. □ 

In the following lemmas we will pass without comment between talking about 
terms or sets of terms in submatrices of a matrix, and subgraphs or cycles 
in the corresponding SR graph. The next formula gives us a way of deciding 
whether two terms in a determinant expansion have the same sign. 

Lemma 5 Consider any square submatrix S{'y\6) of a matrix S. Consider 
any two nonzero terms and Tp in the determinant expansion of S[^\S\ 
corresponding to permutations a and (5 of 5. Then 




(4) 



where \Ce\ is the number of e-cycles in EaU Ep. 



PROOF. Let k 



— I7I — \S\. By definition 



k 



T^Tp^P{a)P{P)l[S,^^^S,^p^. 



i=l 



Let 9 be the set of indices for which ccj 7^ We can write 



ie{i,...,fc}\6> iee 
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When ^ e {1, . . . , k}\e, S.,o..S,^p^ = Sl^^ > 0. So 



sign(TaT^) = P{a)P{(3) sign | H •S'7>a.5'^,ft 



Write j3 o as a product of disjoint nontrivial cycles, and identify each of 
these with a cycle in G(7|5). Let the set of o-cycles in this set be Co and the 
set of e-cycles be Ce, with C = CoUCe- Associate with each cycle c e CoUCe the 
corresponding index set c, i.e., i E c 4^ S^^a,, S^^p. G c. Thus corresponding to 
the sets Co and Ce are the sets of index sets Co and Ce- Since any two cycles 
are edge-disjoint, Co U Ce is a partition of 6, and we can define 



^„=Uc, 9^=1)^ with \eo\ = Y.\cl l^el = E 

ceCo CGCe CGCo CGCe 

Clearly 9 = 9oU 9e. We can write 



c 



^CGCo«ec / \cGCe*ec 



So 



sign(T„T^) =P(a)P(/5) [] sign n^7,«,^7.ft 11 sign [] ^7, -.^7,^ 

=p(a)p(/5)(n(-i)i^|-M f ii(-i)'^' 

\cGCo / \cGCe 

= P(a)P(/5)(-l)l''°l+l'''=l-l'^°l 
= P(a)P(/3)(-l)l''l-l^°l. 

Applying Lemma H] to (3 o gives us that 

P{a)P{(3) = P{(3 o a-^) = (-l)l^l-l^l , 

so that 

sign(T„T^) = = (_i)2|e|-|c|-|c„| ^ (_iyc|+|c.| ^ (_^yc.| 

This proves the result. □ 
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A corollary of the previous lemma is that when the intersection of two term 
subgraphs contains only o-cycles then the two corresponding terms have the 
same sign. 

Lemma 6 Consider any square suhmatrix >S'(7|5) of a matrix S. Let a and (3 
he permutations of 5 such thatT^ andTp are nonzero terms in the determinant 
expansion of S['y\S]. If all cycles in U are o-cycles, then sign(Ta) = 
sign(r^). 



PROOF. Since in this case, there are no e-cycles in E^UE^, applying Eq. (jlj) 
immediately gives 

sign(T„r^) = (-l)l'^^l = (-l)° = l. 
Thus sign(TQ,) = sign(T/3). □ 



Example. Consider the matrix and corresponding SR graph shown in Fig- 
ure [21 Observation of the SR graph tells us that the determinant expansion 
of the matrix contains two terms, and these are of the same sign. Details are 
provided in this figure legend. 
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b 





M = 


— c 





d 







— e 
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Fig. 2. Left. A matrix M. Assume that a, b, c, d,e, f > 0. Right. The corresponding 
SR graph. Identifying edges in the graph with their labels, the graph contains two 
term subgraphs, Ei = {a,d,e} and E2 = {6, c, /}. As Ei U E2 contains a single 
o-cycle [a, b, e, f, d, c], we can infer from Lemma [6] that the corresponding terms Ti 
and T2 in the expansion of det(M) satisfy sign(Ti) = sign(r2). Indeed it is easy to 
confirm that det(M) = Ti +T2 where Ti = ade and T2 = bcf. 



The next lemma shows that having cycles which are both e-cycles and s-cycles 
in a graph means that some terms in a determinant expansion sum to zero. 

Lemma 7 Consider any square suhmatrix 5'(7|(5) of a matrix S. Let a and (3 
he permutations of 5 such thatT^ andTf^ are nonzero terms in the determinant 
expansion of S[^\S\. Assume that E^ 1-) E^ contains exactly one cycle C, and 
this cycle is both an e-cycle and an s-cycle. Then + Tp = 
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PROOF. Let I7I = \5\ = k. By definition 



k k 

T^+Tp = p{a) n + pm n ^7,^ 

i=l i=l 

As usual, let 9 be tlie set of indices for wfiicli a, 7^ I3i so that by assumption, 
5*^,0,, S^.ft e C iff i e 6'. Defining Ci = {S^^a,}iee and C2 = {^^.ftjiee gives us 
a disconnecting partition of C. We can write 



+ = n ^7.«. ^(«) n 57,<.. + p{f3) n -57.^ 

\je{i,...,fc}\e / \ iee iee 

= p{a) I n ^7.«, 1 1 n + o a-') n ^^.^ 

\ie{i,...,A;}\6i / \iee iee 



= P{a) n S,^^^ (sign(Ci)val(Ci) + P(/3 o a-i)sign(C2)val(C2; 
Vje{i,...,fc}\e / 

(3 o can be written as a single cycle of length \9\, and so from Lemma HI 
P(/?oa-i) = (-1)1^1-1. Le., 

T„+T^ = P(«) I n "S"-,.,. I (sign(Ci)val(Ci) + (-l)l^l-isign(C2)val(C2)) 

Since C is an e-cycle we have 

sign(C2)/sign(C0 = sign(Ci)sign(C2) = sign(C) = (-1)1^1. 
Substituting into the expression for + T^, we get: 

T^ + Tp = P{a) I n S,^-, I sign(Ci) (val(Ci) - val(C2)) . 
\ie{i,...,k}\(> J 

However since C is an s-cycle, val(Ci) — val(C2) = 0, giving Tq, + T^j = 0. □ 



Example. Consider the matrix and corresponding SR graph shown in Fig- 
ure [3l Observation of the SR graph tells us that the matrix is singular. Details 
are in the figure legend. 

The next lemma is a consequence of the fact that if a matrix is nonsingular, 
then not all terms can be paired off in the way carried out in Lemma [71 

Lemma 8 Consider any nonsingular square submatrix S{'~f\6) of a matrix S . 
Let a and [3 he permutations of 5 such that T„ and Tp are nonzero terms in 
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M 



—a 


h 





—c 





h 





— c 


a 




Fig. 3. Left. A matrix M. Assume that a,b,c > 0. Right. The corresponding 
SR graph. There are two term subgraphs Ei = {S1 — R2, S3 — R3, S2 — R1} and 
E2 = {Si — Ri, S'a — i?2, 5'2 — i?3}. As EiU E2 contains a single e-cycle which can 
be computed to be an s-cycle, Lemma [7] imphes that that the corresponding terms 
Ti and T2 in the expansion of det(M) satisfy Ti +T2 = 0, and hence that M is 
singular. Indeed it is easy to confirm that det(M) = Ti + T2 where Ti = abc and 
T2 = —abc. 

the determinant expansion of S [-fid]. Assume that U Ep contains exactly 
one cycle C, which is both an e-cycle and an s-cycle. Define Ci = C (1 Ea and 
C2 = C n Ei3 so that {Ci,C2} is a disconnecting partition of C. Then S['y\6] 
must contain a term T„ such that Ci (f_ E„ and C2 E„. 



PROOF. Assume the contrary. Take any term To- in S'[7|5]. If E^ contains all 
the edges from Ci, then we can construct a new term subgraph Er = (-Eo-\Ci)U 
C2 and a corresponding term Tj. in 5'[7|(5]. Alternatively if E„ contains all the 
edges from C2, define E^. = {E„\C2) U Ci with corresponding term T^. By 
construction, E^r U Er contains only one cycle which is an e-cycle and an s- 
cycle and so, by Lemma [3, T„-\-Tr = 0. Thus all terms in 5'[7|(5] fall into pairs 
which sum to zero and 5(715) is singular. □ 



The next lemma tells us a fact which is geometrically obvious about how a 
term subgraph can intersect a cycle: either it contains all members in one half 
of a disconnecting partition of the cycle, or it contains an edge not in the cycle 
but incident on a vertex in the cycle. 

Lemma 9 Consider a square matrix S and the associated SR graph G. As- 
sume that G contains a cycle C. Let {Ci,C2} be a disconnecting partition of 
C , and E be any term subgraph in G. LfCi E and C2 <^ E, then E contains 
an edge incident on an S-vertex in C but not itself in C . 



PROOF. Let |C| = 2k. Since ii^ is a term subgraph in G, it has an edge 
incident on each vertex in C . Either some of these are not in C, or E contains 
k edges in C. In the latter case, either they are all in Ci, or they are all in 
C2, or two of them share a vertex. The first two possibilities are ruled out 
by assumption, and since is a term subgraph, the third situation is not 
possible. So E must contain an edge e incident on a vertex in C, but not itself 
lying in C. If this is an S-vertex we are done. If it is an R- vertex, but e is also 
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incident on another vertex in C, then this second vertex must be an S- vertex, 
and again we are done. Otherwise e is incident on an R-vertex in C and some 
vertex outside C, in which case E\{e} contains edges incident on k S-vertices 
and k — 1 R- vertices in C . There must hence be some edge e' e E\{e} incident 
on an S- vertex in C but not itself in C. □ 



Note that in the statement and proof of the above lemma we could exchange 
"S- vertex" and "R-vertex". 

Our final preliminary lemma tells us that if an S-to-R path "slices" an e-cycle 
in a particular way, this implies the existence of two e-cycles with S-to-R 
intersection. 

Lemma 10 Consider a matrix S and the associated SR graph G. Assume 
that G contains an e-cycle C . Assume that there is a path D in G joining an 
S-vertex in G to an R-vertex in G, but such that D and G are edge- disjoint. 
Then G contains two e-cycles with S-to-R intersection. 



PROOF. Let G be of length 2r with G = {Saibi, Sa^bi, Sa2b2, Sa^br, ^'a^fej, 
and let D join S-vertex aj to R-vertex bk- Decompose G into the two edge- 
disjoint paths G = GiU G2 where 

Gl = {Sajbj, Saj_f^-^bj, ■ ■ ■ , Sa^bk}: ^2 = {Saf.^ibky Saf._^.j^bk+iy • • • ) Sajbj^i} ■ 

Note that |Ci| and IC2I are both odd since Gi and G2 are S-to-R paths. The 
situation is illustrated in Figure HI 



Fig. 4. A schematic representation of the situation in Figure [TOj Ci and C2 are 
S-to-R paths which partition C. D joins an S-vertex to an R-vertex in C but is 
edge-disjoint from C. (Note that in this figure, bold lines indicate edges of arbitrary 
sign, while dashed Unes indicate paths of arbitrary length.) 



There are two cycles Di = D U Gi and D2 = D U G2 such that Di n D2 = D, 
Di O G = Gl and H C = G2. We will show that Di and D2 have opposite 
parity and hence one of them must be an e-cycle. Note that |Ci| -|- IC2I = 2r, 
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\Di \ = \D\ + \Ci\ and ID2I = \D\ + IC2I. Subtracting the last two expressions 
gives ID2I - l^il = IC2I - \Ci\ = 2(r - |Ci|). 

The fact that C is an e-cycle means that either 

(1) r is odd and sign(C) = — 1, or 

(2) r is even and sign(C) = +1. 

Case 1. Since sign(C) = sign(Ci) sign(C2) = —1: 

sign(Ci) = -sign(C2) 

Now r and |Ci| are odd, so r — |Ci| is even, so I-D2I — \Di\ = 2(r — |Ci|) is a 
multiple of 4, i.e., (-l)l^2|/2 _ (_i)Pi|/2_ xhen 

P{Di) = (-l)l^^l/2sign(Ci)sign(D) = (-l)l^^l/2[-sign(C2)]sign(Z}) = -P(Z^2). 

So P{D2) = —P{Di) and one of Di or D2 must be an e-cycle. 

Case 2. This time sign(C) = sign(Ci) sign(C2) = 1 so 

sign(Ci) = sign(C2) 

Now r is even and \Ci\ is odd, so r — |Ci| is odd. As a result ID2I — |-Di| = 
2(r - \Ci\) is not a multiple of 4, and (-1)1^21/2 = (_i)Pi|/2+i_ Again we get 
P{D2) = —P{Di) and one of Di or D2 must be an e-cycle. 

In each case, one of Di or D2 is an e-cycle. Moreover both Di and D2 intersect 
C along an S-to-R path (either Ci or C2). As this is the unique component of 
their intersection, they have S-to-R intersection. □ 



5 Relationship between sign nonsingularity and o-cycles 

There is a very simple and elegant relationship between properties of square 
submatrices of a matrix and the non-existence of e-cycles in the corresponding 
SR graph. The results we prove in this section are weaker than our main result 
in the next section, but have a certain generality to them. At an abstract level, 
the results simply describe an interesting relationship between a matrix and 
an associated bipartite graph. Comments in Chapter 3 of [4j suggest that 
these results may be known or suspected, but we have not found a proof in 
the literature. In any case they are an easy corollary of the previous lemmas. 
First, analogous to the definition of sign- nonsingularity, we define a square 
matrix S to be sign singular, if all matrices with the same sign pattern as 
S are singular. 
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Theorem 11 Consider a matrix S and the associated SR graph G. If all 
cycles in G are o-cycles, then all square submatrices of S are either sign non- 
singular or sign singular. 



PROOF. Consider any square submatrix S{'~f\6) of 5* with rows and columns 
indexed in the usual way by ordered sets 7 and 6 respectively, where I7I = 
1^1 = k. Consider permutations of 6 and the corresponding terms in S'[7|5]: 

1) If all such terms are zero then the S{'y\S) is sign singular. 

2) If there is a single a for which is nonzero then 5'[7|(5] = so clearly 
S{j\S) is sign nonsingular. 

3) Consider any pair of nonzero terms in >S'[7|5] corresponding to permutations 
a and (3 of 6. Since all cycles are odd, Lemma [6] gives us that sign(Ta) = 
sign(T/3). Since a and (3 were arbitrary, this means that all nonzero terms in 
5* [7 1 5] have the same sign. □ 



The previous theorem has a converse: 

Theorem 12 Consider a matrix S and the associated SR graph G. If all 
square submatrices of S are either sign nonsingular or sign singular, then all 
cycles in G are o-cycles. 



PROOF. Assume that G has an e-cycle C of length 2r including S-vertices 
from a set 7 = {71, 72, ... , 7^} and R-vertices from a set 5 = {61, 62, ■ ■ ■ , Sr}. 
There is some permutation a of 5 such that C consists of the edges 



U {s- 



J- 



iG{l,...,r} 



Setting Pi = ctj+i defines a permutation /? of 6. Clearly S{'y\6) is not sign 
singular since 

T« = P(a) n and = P(/3) fj 5,,^, 

i=l i=l 

are nonzero terms in 5'[7|(5]. Since E/^ = C, and C is an e-cycle, applying 
Eq. (Hj) gives: 

sign(T„T^) = (-l)i = -l. 
As Tq, and Tp have opposite signs, 5'(7|(5) fails to be sign nonsingular. □ 
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Note that in the terminology of [3] a matrix which is either sign nonsingular 
or sign singular is a matrix with signed determinant. We can state the 
previous two theorems as a single result: 

Corollary 13 Consider a matrix S and the associated SR graph G. The fol- 
lowing two statements are equivalent: 

(1) All square submatrices of S have signed determinant. 

(2) All cycles in G are o-cycles. 



PROOF. This follows immediately from Theorems [TT] and [T^ □ 



Remark. Corollary [13] is a general statement about matrices and the corre- 
sponding bipartite graphs. Note that insisting that all square submatrices of a 
matrix have signed determinant is considerably more restrictive than insisting 
that they are all either sign nonsingular or singular (i.e. that the matrix is 
SSD). This is because it is a frequent occurrence for a square matrix to be 
singular while failing to be sign singular. 



6 A graph-theoretic condition ensuring injectivity 

Define the following condition on any SR graph: 

Condition (*): all e-cycles in the SR graph are s-cycles, and no two e-cycles 
have S-to-R intersection. 

In [7j it was shown that Condition (*) on the SR graph of a reaction system 
with mass-action kinetics and with outflows is sufficient to ensure injectivity 
of the system. For NIC reaction systems this means, by results in [3], that 
Condition (*) is sufficient to ensure that the stoichiometric matrix S of the 
system is "weakly sign determined" (WSD), a less restrictive condition on S 
than the requirement that it should be SSD. However results in ^ and [3J left 
open the question of whether Condition (*) implies that the stoichiometric 
matrix is SSD giving injectivity for essentially arbitrary kinetics. Below we 
show that this is indeed the case - Condition (*) implies that the stoichio- 
metric matrix is SSD. Since SSD implies WSD, as a corollary we reproduce 
the result that can be inferred from ^ and [3J . We also show by example that 
Condition (*) is not a necessary condition for a matrix to be SSD. 
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6.1 Condition (*) ensures that a matrix is SSD 

Theorem [1] is the main result of this paper: that an NIC reaction system 
whose SR graph fulfils Condition (*) has stoichiometric matrix which is SSD, 
and is hence, with the outflow conditions detailed at the beginning, incapable 
of multiple equilibria. 

PROOF of Theorem [1]. Assume that S is not SSD, i.e., there exists some 
square submatrix 5'(7|5) of S which is neither sign nonsingular nor singular. 
Since 1x1 matrices are trivially sign nonsingular or singular, I7I = \5\ > 2. 
The proof will proceed by showing that the corresponding subgraph G{'~f\6) 
either contains an e-cycle which fails to be an s-cycle or contains two e-cycles 
which have S-to-R intersection. 

Firstly, it is immediate from Theorem [11] that 0(^16) contains an e-cycle. 
If I7I = 2, then there are exactly two terms in 5'[7|5], i.e., S'^^^^S'^j^j and 
—5^152 5*^2 <5i- Since S{'y\S) is not sign nonsingular, 

sign (5*^^ 5^5*^2 52) = sign(S'^i52'S'72<5i) • 

Since 5'(7|5) is not singular, 5'[7|(5] = S'^^5jS'^2<52 ~ 'S'^i<52'S'72<5i 7^ 0. Hence, 
by the definition of an s-cycle, the e-cycle in G{'y\6) consisting of the edges 
{S^j^Si, 5'^i(52 5 S-y^s2: 'S'^2<5i} is not an s-cycle and Condition (*) is violated. 

So now assume I7I > 3. Consider two terms Tq, and of opposite sign in 5'[7|(5] 
corresponding to subgraphs and Ejj in G. By LemmaEl EaUEyj contains an 
e-cycle. Take an e-cycle C C {EaUE.ri) and construct a disconnecting partition 
of C into Ci = CnEa and = CnE^ according to Eq. Let |C| = 2k for 
some 2 < /c < I7I so that \Ci\ = IC2I = k. Define Ejs = {Eo,\Ci) U C2. Ejs has 
been constructed so E^ U Ep contains only one cycle (i.e., G). Corresponding 
to Ea and Ep are nonzero terms T„ and in S'[7|5]. Bearing in mind that 
there is only one cycle in Ea U Ej3 and this is even, Eq. (jl]) gives: 

sign(T„T^) = (-1)1 = -1. 

There are now two cases to consider. If + 7^ 0, then vaA{Ea) 7^ val(i?/3). 
I.e., val((E„\C) U (E„ n C)) ^ val((E^\C) U{E^nG)). Since E„\C = E^\C, 
we get val{Ea fl C) 7^ val(-E^ fl G). Since EaH G and E13 n G together make 
up a disconnecting partition of G, this means that G is not an s-cycle and we 
are done. 

So assume that Ta + T13 = 0. Now by Lemma [8] we can find a permutation a 
of 6, and corresponding term T„ in 5'[7|(5] such that Gi (f. E„ and G2 ^ E„. 
As a result, by Lemma [HI E^^ contains some edge incident on an S- vertex in 
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C, but not itself in C. In other words, there is an index q such that edges 
S-rqaq, S^^i3g e C, but aq ^ aq and a, 7^ Pq so that S^^^^ ^ C. 

Consider a o as the product of disjoint cycles. Since aq 7^ aq, a nontrivial 
cycle involves aq. Further, since aq 7^ Pg, the corresponding cycle in G is 
distinct from C but intersects C at S-vertex 7^. Let this cycle be termed Caa 
(to remind us that it is composed of edges from Ea and E„) and have length 
2r (2 < r < I7I). In the usual way, we can follow Caa- there are distinct indices 
ai, . . . , a,,. G 7 and bi, . . . ,br G 5 such that the edges Sa b occur in Ea and 
occur in E^j. Assume (without loss of generality, i.e., by reordering the 
sets tti and hi if necessary) that ai = 7^ and 1)2 = a q, so that Sa^h2 is the edge 
incident on an S-vertex in C, but not itself lying in C. 

Follow the cycle Caa starting at Sa^ij2 ^ T„ 

Sa\h2 ^ -^<T) S(i2h2 ^ -^ct) ^a2b3 ^ -^<t) ■ ■ ■ 

Since this is a cycle, eventually some vertex from the sequence of alternating 
R- and S- vertices indexed by (62, ^2, 63, . . .) must be a vertex in C. But this 
cannot first happen at an S-vertex. Suppose the contrary and S-vertex aj is in 
C, while R-vertex bj is not. Since edge Sajbj is in Tq,, but not in C, this implies 
that S-vertex aj has three edges from Ta U incident on it - the two edges 
in C along with the edge Sa b - But this is impossible from the discussion 
in Section [31 This means that from the vertex sequence (62, ^3? • • •) the 
first vertex to lie in C must be an R-vertex. Let this vertex be bj {j = 2 is 
possible). Define the path D = {Saj^b2, Sa2b2, ^'aabg, • • • , 5'a^_i6j}- D is an S-to-R 
path starting and terminating at vertices in C but edge-disjoint from C. 

Now applying Lemma [10] shows that there are two e-cycles in G (one of which 
is G and one of which is made up of D and some part of G) which have S- 
to-R intersection. Thus if the stoichiometric matrix S is not SSD, then the 
associated graph G necessarily fails Condition (*). The result is proved. □ 



6.2 Condition (*) is not necessary for SSD (or WSD) 

We present an example to illustrate that Condition (*) on the SR graph is not 
necessary to give a system with an SSD matrix. Consider the system of three 
reactions 



D ^ A + B + C, E ^ A + B + C, F^A + B. (5) 
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This system has stoichiometric matrix 



S 



111 
111 
1 1 
-10 
0-1 
0-1 



which can quickly be computed to be SSD, implying that the system with 
inflows and outflows forbids multiple equilibria. The SR graph for the system 
is shown in Figure [51 Although all e-cycles are s-cycles, there are e-cycles with 
S-to-R intersection, and thus the graph fails Condition (*). 




Fig. 5. The SR grapti of reaction system ([5]). Edge labels are all 1 and have been 
omitted. As a consequence all cycles are s-cycles. However, there are a number 
of e-cycles with S-to-R intersection, for example the cycles A — R2 — B — and 
A — R2 — C — R1 intersect along the S-to-R path A — R2. 

Thus Condition (*) is not necessary for injectivity in the case of a general sys- 
tem of NIC reactions. Since the condition that S is SSD is stronger than the 
condition that S is WSD, clearly Condition (*) is not necessary for injectivity 
in the case of a system of mass action reactions. 



7 Conclusions 



We have described several new results for the class of NIC reaction systems 
with arbitrary chemical kinetics. The key has been to associate with any real 
matrix a signed, labelled, bipartite graph termed an SR graph. Given any 
matrix S and associated graph G we have shown that: 

(1) All cycles in G are o-cycles iff all square submatrices of S have signed 
determinant (Corollary [T3D . 
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(2) Condition (*) on G is sufficient, but not necessary, to guarantee that 5* 
is SSD (Theorem [H and the counterexample in Section [6^ . 

If S is the stoichiometric matrix of a chemical reaction network, and we assume 
the outflow conditions in System ([2]), either graph-theoretic condition implies 
immediately that multiple equilibria cannot exist [3]. Otherwise, they imply 
that multiple positive nondegenerate equilibria cannot exist [6]. These results 
apply to large classes of realistic biochemical reaction networks; for example, 
they apply to the enzymatic reaction networks discussed in [8], without the 
assumption of mass-action kinetics. There are natural further extensions of 
this work to the situation where the NIC condition is dispensed with. These 
topics are pursued in |2j. 

The identification of a condition on G which is equivalent to S being SSD 
remains an open - and probably very difficult - problem. However, there are 
certain ways in which G can fail Condition (*) which ensure that S fails to be 
SSD. These will be discussed in future work. 

Finally, it should be mentioned that there are interesting relationships between 
the results in this paper and approaches to injectivity involving so-called "in- 
teraction graphs". The approach taken in [15] is most similar to that taken 
here, with the key difference that interaction graphs and SR graphs are for- 
mally quite different objects. The details of how interaction graph approaches 
and SR graph approaches are related to each other will be explored in future 
work. 
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